Probing the quantum nature of spacetime by diffusion 
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Many approaches to quantum gravity have resorted to diffusion processes to characterize the spec- 
tral properties of the resuhing quantum spacetimes. We critically discuss these quantum-improved 
diffusion equations and point out that a crucial property, namely positivity of their solutions, is 
not preserved automatically. We then construct a novel set of diffusion equations with positive 
semi-definite probability densities, applicable to Asymptotically Safe gravity, Hoi^ava-Lifshitz grav- 
ity and Loop Quantum Gravity. These recover all previous results on the spectral dimension and 
shed further light on the structure of the quantum spacetimes by assessing the underlying stochas- 
tic processes. Pointing out that manifestly different diffusion processes lead to the same spectral 
dimension, we propose the probability distribution of the diffusion process as a refined probe of 
quantum spacetime. 

PACS numbers: 02.50.-r, 04.60.-m, 05.40.-a, 05.45.Df, 05.60.-k 



I. INTRODUCTION 

In quantum gravity it is commonly expected that quan- 
tum, fluctuations will lead to a spacetime structure that 
departs from a smooth classical manifold at short scales. 
A common theme shared by many quantum gravity can- 
didates focuses on the characterization of the result- 
ing quantum spacetimes and possible observable conse- 
quences. In this context, the analysis of the spectral 
properties, which can be probed by the diffusion of a fidu- 
cial test particle on the effective quantum spacetime, have 
turned out to be quite useful. In particular, they allow 
discriminating between the pictures of quantum space- 
time emerging from different theories. 

The simplest spectral quantity carrying nontrivial in- 
formation of the spacetime structure is the spectral di- 
mension ds ■ This notion of dimensionality, which is dis- 
tinct from the topological or HausdorfF dimension, arises 
in the context of a diffusion process: From the prob- 
ability density P{x,x',a) of a diffusing particle on a 
background, one can define a return probability V{a) = 
V^^ J P(x,x,a). Herein, x and x' denote coordinates 
on the Euclidean spacetime, V is the volume and a an 
external diffusion time. The spectral dimension for a 
background with fixed dimensionality is then defined aqj 



da = — 21im 



d\xiV{a) 
dlna 



(1) 



By now, the spectral dimension has been determined 
within several quantum gravity approaches, including 
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^ At this stage the derivative symbols are not strictly necessary 
due to the presence of the limit, but will become important in 
the more general case of a scale-dependent spectral dimension. 



Causal Dynamical Triangulations (CDT) in four [l|, [2| 
and three spacetime dimensions [3, |j|j Euclidean Dy- 
namical Triangulations (EDT) [5| , Loop Quantum Grav- 
ity (LQG) 0], Asymptotically Safe gravity (or Quantum 
Einstein Gravity, QEG) \M\^ , nonlocal gravity [11| , and 
also Hofava-Lifshitz (HL) gravity [IJ, [ij]- Remarkably, 
many of these works assess that ds = 2 at microscopic 
scales |14| - [la |. Subsequently, it actually turned out to 
be useful to generalize the definition ([1]) by dropping the 
limit cr — > and considering ds{cr) as a function of the 
diffusion time. This generalization effectively allows one 
to characterize a spacetime with multifractal structures 
where the spectral properties change on various typical 
length scales, e.g., when considering the spacetime struc- 
ture at classical, Planckian, and sub-Planckian distances. 

To arrive at the return probability, two different routes 
have been followed in the literature. The Monte Carlo ap- 
proaches to quantum gravity, foremost CDT and EDT, 
approximate the quantum spacetime by a piecewise lin- 
ear manifold built from d-dimensional simplices [ITj . In 
this setting, the quantum nature of spacetime is encoded 
in the gluing of the building blocks. The return proba- 
bility is then obtained by studying a standard random 
walk on the resulting discrete quantum spacetime. In 
other words, a classical random walk is used to capture 
the fractal features of the geometry. The probe parti- 
cle thereby has equal probabilities for moving from its 
current simplex to one of its neighbors in each time step 
[ll y-Q- This construction implies that, after each time 
step, the probe particle has to be located somewhere on 
the geometry, guaranteeing a positive semi-definite prob- 
ability density P{x, x' , a). 

Within analytical approaches to quantum gravity, it 
has been proposed that quantum effects in spacetime, 
foremost a dynamical dimensional change, can be cap- 
tured by a modification of the Laplacian operator ap- 
pearing in the classical diffusion equation |6l-[l3l [l8| . We 



will call this new operator "generalized Laplacian" or, 
as done in probability theory, spatial generator. Thus, 
one uses an anomalous diffusion equation on classical flat 
spacetinic to mimic the quantum structure. The under- 
lying physical picture is that the effective metric "seen" 
by the diffusing particle actually depends on the momen- 
tum of the probeo Expressing this metric through a fixed 
reference scale leads to a modified diffusion equation pro- 
viding an effective description of the propagation of the 
probe particle on the quantum gravity background. For 
Asymptotically Safe gravity jl9l - l22l | in d dimensions, this 
procedure gives rise to a multifractal structure where the 
spectral dimension is of the form [9| 



2d 



(2) 



Here, the parameter S captures the quantum effects and 
actually depends on the probed length scale. One can 
identify three characteristic regimes where the spectral 
dimension is approximately constant over many orders 
of magnitude. At large distances, one encounters the 
classical regime where 6 = and the spectral dimension 
agrees with both the Hausdorff and topological dimension 
of the spacetime. At smaller distances one first encoun- 
ters a semiclassical regime with 5 = d, before entering 
into the fixed-point regime with 5 = 2. 

Let us clarify some conceptual points arising in this 
approach of characterizing the quantum spacetime. Cru- 
cially, diffusion processes probe the properties of the ef- 
fective background spacetime only. There is no relations 
to fluctuations of the backgrounds which would be en- 
coded in the graviton propagator. Accordingly, in our 
context the spectral dimension is probably not directly 
understood as "the dimension of momentum space," as 
conjectured in fractal geometry [23, ^M- Furthermore, 
the diffusing particle is a fictitious probe, and no back- 
reaction of the particle on the background is included. 
Finally, the diffusion time a is an external time, thus it 
is not related to a (causal) propagation of a matter field 
on the quantum spacetime. In other words, the scaling 
of the graviton and also matter propagators in quantum 
gravity need not show an effective two-dimensionality, 
even in the case ds = 2 |25l - [2g | . We thus stress that the 
spectral dimension is a useful tool to encode properties 
of an effective spacetime in quantum gravity, but it is 
not necessarily related to the propagation of a physical 
particle. 

A potential issue related to the use of the modified 
diffusion equations proposed in the quantum gravity lit- 
erature [61, [7|, i^, ilCi, il2|, ilSj] is their potential risk of giving 



^ Although a scale-dependent metric does not seem to make sense 
at a first glance, one should note that a unique, dynamically 
determined metric emerges as the solution to the full quantum 
equations of motion in the infrared. Determining momentum 
scales with respect to this metric then gives meaning to the no- 
tion of a scale-dependent family of metrics. 



rise to "negative probabilities." In fact, we will explic- 
itly demonstrate in the next section that the solutions of 
these modified equations are not positive semi-definite, 
so that the interpretation as probabilistic processes is 
lost. We take this observation as a clue that these de- 
scriptions are incomplete, and formulate several propos- 
als for restoring the probabilistic nature of the modified 
diffusion processes in Sees. IIIII to IVII This allows us to 
make another step forward in the exploration of quan- 
tum spacetimes: The diffusion probability encodes more 
information on the underlying spacetime than the spec- 
tral dimension alone, and permits to study the quantum 
properties of spacetime in more detail. In particular, 
we provide a setup which allows one to construct diffu- 
sion equations accommodating generic nontrivial scaling 
properties even within an anisotropy between space and 
time. We thus expect our setup to be relevant to several 
quantum gravity approaches, encompassing Asymptotic 
Safety, Hof ava-Lifshitz gravity and Loop Quantum Grav- 
ity, as well as the tentative continuum limit of CDT and 
EDT. 



II. DIFFUSION EQUATIONS WITH NO 
PROBABILISTIC INTERPRETATION 

Non-relativistic Brownian motion in d dimensions with 
initial condition P{x, x' , 0) = 6{x — x') is described by 



{d^~Vl)P{x,x',a) =0. 



(3) 



One way to encode quantum-gravity effects is to modify 
the Laplacian V^ ~ .g^'^V^Vjy by replacing the classical 
spacetime metric g^^^ by a suitable quantum object [7|, 



(a,-(A))P(.T,.T',a)-0, 



(4) 



where (A) is a generalized Laplacian. In Asymptotically 
Safe gravity, this object is naturally constructed from the 
vacuum expectation value of the quantum metric which 
can be computed from the scale-dependent flowing ac- 
tion Ffc. Similarly, heuristic arguments in Loop Quantum 
Gravity allow one to obtain a quantum version of g^i^ 
based on the area operator [6| , while anisotropic Hof ava- 
Lifshitz gravity resorts to a Laplace operator that con- 
tains higher powers of spatial derivatives [l^, ^M, ■ 

The solution P for these cases can be evaluated 
by resorting to a flat reference spacetime at macro- 
scopic scales. One can then express P(x,a;',cr) through 
its Fourier transform. For example, the application 
of Renormalization Group (RG) improvement schemes 
within Asymptotic Safety yields a solution of the form 
0i 



d'^p 



P,,o(x,x',a)= / ,^4Le^P-(--')e-P^^(^^) 



(27r) 



(5) 



The function F{p^) = Ap/A^^, encoding the quantum 
behavior, is determined by the value of the dimensionful 



cosmological constant seen at energy scale 



' relative 

to its infrared (IR) value. The short-distance behavior 
found in the LQG-inspired calculation of [y] is a special 
case of (O with Fjp^) = p^ . 

In HL gravity [30| , the quantum effects are encoded in 
a modified dispersion relation uP = f{p^)- Herein the 
dependence of the function / on the spatial momenta p 
is parameterized by 



f{f)=P' 



l + Bf + Cp^ 



(6) 



The parameters B,C, D have been determined through 
a comparison to three-dimensional CDT data Q and in 
the following we will use the "preferred values" B = 
-1.18, C = 344.47,1? = 10.08 found in jjj. The 
anisotropic dispersion relation then results in the expres- 
sion 



PnLix,x\a) 



(2^)3 



(7) 

Here the Wick-rotated time coordinate t of the manifold 
is not to be confused with the external diffusion time a. A 
rather generic consequence of having nontrivial functions 
F{p^) and f{p^) is the lack of an a priori guarantee of 
a positive semi-definite solution: in general, the Fourier- 

2 + 5 

transform of a non-Gaussian function of the form e""'^ 
is not positive semi-definite. For illustrative purposes, we 
exemplify this feature based on the expressions ([S|) and 
dni). Since the functions F and / depend on the abso- 
lute value of the Euclidean momentum only, the angular 
integrations in momentum space can be carried out by 
applying 



d'pf{p)e 



{2n) 



dpp^Ja_^{pr)f{p), (8) 



where Ji,{x) is the Bessel function of the first kind and 
order v and r = \x\. In this way, the probability densities 
reduce to the Hankel transform of the modified exponen- 
tials. Explicitly, 



PQEG{r,a) = 
while 



(2^'-) 2 Jg 



dpp'^'Jd/2-i(pr)e-''P'^<^P'\ 



(47rcr)i/2 7g 2tt 



kdk 



(9) 



Jo(fcr)e-"-''(^'). (10) 



The remaining integral in Eq. (jlOp can be evaluated 
numerically. For fixed values of a and t, snapshots of 
Eqs. (O and ([10]) arc shown in Figs. [H [2l and IH As 
illustrated by these plots, a common feature shared by 
all these modified diffusion equations is that the result- 
ing probability densities for the probe particle are no 
longer positive semi-definite. Thus, strictly speaking, 
the "quantum-improved" diffusion equations do not give 




FIG. 1. The function PqEcCr, cr = 1), Eq. ((9]), evaluated for a 
typical RG-trajectory of Quantum Einstein Gravity in d = 3. 
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FIG. 2. The function PqEcCr, a = 1), Eq. ((gj, evaluated for 
the ultraviolet (UV) limit in LQG-inspired scenario of [y|, by 
setting d = 4 and F{p^) = p^. 
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FIG. 3. The function PHL(r,t = l,o- = 1), Eq. ^, obtained 
from the anisotropic diffusion equation of HL gravity [l3l |. 
The occurrence of negative values is a generic feature and 
does not rely on the special choice of parameters in Eq. (O . 



rise to a well-defined diffusion process. Accordingly, one 
should proceed with great care when interpreting the re- 
sulting quantity ds as a spectral dimension. Moreover, 
at a conceptual level this feature constitutes a mismatch 
between the spectral properties derived within Monte- 
Carlo studies and the improved diffusion equations. We 
take this as a motivation to investigate the possibility of 
"quantum-improved" diffusion equations which give rise 
to a manifestly positive-definite probability density in the 
next sections. As an additional bonus, the construction 
of these will also allow us to gain further insight into 
the nature of the diffusion process and the underlying 



stochastic interpretation in quantum gravity. 



III. GENUINE DIFFUSION ON QUANTUM 
SPACETIME 

The aim of this section is to construct diffusion equations 
which capture the quantum properties of spacetime while 
admitting solutions that are manifestly positive semi- 
definite. Essentially this can be achieved in two ways: 
either by considering diffusion in nonlinear time and/or 
modifying the quantum diffusion equation ^ through a 
nontrivial source term. Both cases will be discussed be- 
low. For concreteness, we will focus on diffusion processes 
in Asymptotic Safety, but the results are also relevant for 
the UV limit of the LQG-inspired model of Ref. [i|. The 
generalization to HL gravity will be discussed in Sec. IIVI 
As a nontrivial finding, our novel diffusion equations re- 
cover all previous results on ds(a). 



The relation between k and a is then fixed on dimensional 
grounds. Since kx is dimensionlcss, Eq. (|12p implies that, 
dimensionally, a ^ k^^^^ and suggests the scale identifi- 
cation 



k ^ a «+2 , 



(13) 



where the proportionality constant has been absorbed 
into diffusion time a. By changing the diffusion time- 
variable from a to a^ with 



(14) 



(5 + 2' 



the resulting equation can be cast into a diffusion equa- 
tion in nonlinear time: 



(^^'^'l^^^'"''") 



(15) 



The probability density resulting from this diffusion 
equation is given by a Gaussian in r = jx — x'j: 



A. Diffusion in nonlinear time 

The crucial step in obtaining Eq. (|4]) via an RG- 
improvement is the identification of the cutoff scale k 
with a suitable physical quantity. Such a procedure al- 
lows one to analyze the effects of quantum gravity in 
a variety of single-scale settings [3l[. Since we have 
seen that, generically, any quantum-gravity modifica- 
tion of the form V^ i— >■ (A) will result in functions 
P{x, x' , a) that have no interpretation as probability den- 
sities, we will explore an alternative RG-improvement 
scheme where the scale k is related to the diffusion time 
a. On general grounds, the identification k — k{a) should 
be such that large diffusion times correspond to the IR 
regime A: — > and short diffusion times to the UV regime 
k — > oo. 

The RG-improved diffusion equation based on non- 
linear time can be derived through a modification of 
the original RG-improvement scheme of Ref. [3|- For 
that purpose, we assume that the scale-dependent so- 
lutions of the equations of motion based on the fiowing 
action Tk give rise to a power-law relation between the 
effective metrics at scale k and the reference scale fcg, 

{gnk « 



k'ig' 



see App. [A] Substituting this rela- 



tion, Eq. (U) becomes 

(9, - fc'(5^''V^V,)fcJ P{x,x',<j) - , (11) 
where k is the RG-scale, fco is the IR reference scale and 



9'' 



is the fixed IR-reference metric which we take to 



be the flat Euclidean metric. Since we want to encode 
the scaling effects in the diffusion time, we multiply the 
equation with k~^ so that the diffusion operator becomes 
a standard second-order Laplacian, 



k-'-^-S/l)P{x,x',a) = 0. (12) 



Pir,a) 



1 



(Anaf^y. 



(16) 



Thus it is manifestly positive semi-definite. Moreover, 
the cutoff identification p^ implies that P[r,a{k)] ex k'^ 
has the correct scaling behavior of a diffusion probability 
in d dimensions, giving further justification to the RG- 
improvement procedurclfl 

The spectral dimension resulting from Eq. d 
dependent of a and reads 



ds 



2d 

2 + 6' 



(17) 



Notably, this expression is identical to the one derived 
from RG-improving the Laplacian of the diffusion equa- 
tion [9|. Thus, the two RG-improved schemes give rise to 
the same profile for the spectral dimension. This implies, 
in particular, that the matching between the spectral di- 
mension obtained from the diffusion in nonlinear time, 
Eq. ([T5|) . and the one measured within CDT is identical 
to the one found in Q- Moreover, the positive semi- 
definite function (J16p has the correct qualitative features 
of a probability density measured within CDT, making a 
future comparison between these more detailed spectral 



^ This diflfusion equation is also relevant to the ease of Loop Quan- 
tum Gravity, where the scaling of the area operator at small 
distances Z <^ Zq suggests to deduce a scaling of the metric in 
the form g^^" (I) ~ l~^9'^'^{lo) 111- Following the steps outlined 
above then allows to arrive at II15II as a description of diffusion 
processes in the LQG-inspired model of [y]. However, caution 
should be exercised in drawing any conclusion regarding the spec- 
tral dimension in Loop Quantum gravity, since computations in 
a bottom-up approach (i.e., by placing a random walker in re- 
alistic LOG graphs) can lead to considerably more complicated 
results |32|. 



properties of the underlying quantum spacetimes mean- 
ingful O 

The averaged squared displacement of the test-particle 
implied by ([T6|) is easily found to be 



V /nonlinear time 



2da^ . 



(18) 



where angular brackets denote the expectation value with 
respect to the associated probability density function, 
(/(a;)) = Jd'^xP{x,x',a)f{x). For (3 = 1, this corre- 
sponds to a standard Wiener process. In the case /3 < 1 
the diffusion time passes slower, and slows down further 
throughout the diffusion process. The diffusion process 
is subdiffusive, i.e., the average displacement of the test- 
particle is less than in the diffusion on an ordinary d- 
dimensional manifold. 

There are actually two possible stochastic processes 
underlying Eg . (|15p . One is scaled Brownian motion 
(SBM) [33t35[, i-e., a Wiener process which takes place 
in "nonlinear time." The second is fractional Brown- 
ian motion (FBM) [3J-[33] , which is a stochastic pro- 
cess with correlated increments, or, in other words, non- 
Markovian. Let us note that, since both stochastic pro- 
cesses lead to the same diffusion equation, we cannot dis- 
tinguish them at this point. They do however differ in 
one crucial property: SBM is actually Markovian, while 
FBM is not. Although we do currently not have enough 
data to rule out a non-Markovian stochastic process in 
either Asymptotically Safe gravity or the tentative con- 
tinuum limit of CDTs, this property could be used to 
distinguish between the processes in the future. Both 
processes and their properties (as well as a third twin, 
not likely to be relevant in the present context) are dis- 
cussed in App.|B] There we recall how diffusion processes 
are classified. This helps to develop a more intuitive un- 
derstanding of the properties of quantum spacetime as 
probed by a diffusing particle. For instance, Eq. (|15p 
corresponds neither to trapped subdiffusion (where the 
test particle spends much time in bound states) nor to 
labyrinthine diffusion, where "obstacles" and "holes" are 
interposed by geometry and topology, as is in accordance 
with the expectation from Asymptotically Safe gravity. 
In particular, this does not correspond to diffusion on 
a fractal in a mathematical senseO The alternative in- 
terpretation is that subdiffusion occurs because of a vis- 
coelastic effect: the test probe is dragged by the complex 



While the probability density P{x,x' ,a) can certainly be mea- 
sured on the quantum spacetimes underlying CDT, no such data 
are currently available. 

Transport on fractals is realized by labyrinthine diffusion |34| : 
subdiffusion is caused by the geometric and topological struc- 
ture of the fractal and, so to speak, the test probe lives in a 
"crowded" environment and meets a number of obstacles and 
dead ends. Diffusion on a fractal is, in general, described by 
a diffusion equation with fractional differential operator d^, see 
Eq. 111911, and an x-dependent diffusion coefficient plus friction 
(see |38| and references therein). 



environment of which it is part. Intuitively, the subd- 
iffusion can be understood as quantum fluctuations of 
spacetime slowing down the propagation of the particle 
as compared to a classical background. It is interesting to 
observe that this property is actually common between a 
variety of quantum gravity approaches, whereas no case 
is known in which quantum effects lead to superdiffusion. 
Let us emphasize again at this point that these findings 
do not allow to conclude that a physical particle propa- 
gating on the effective spacetime is subject to any similar 
effect. 



B. Diffusion employing fractional derivatives 

An interesting alternative to the cutoff identification (J13p 
relates the fc-dependence of Eq. (fT2|) to the order of the 
derivative operator. The resulting diffusion equation is 
formulated in terms of fractional derivatives 



(19) 



yl)Ppix,x',a)=0, 



where d^ is the left Caputo derivative (see, e.g., [39|-|4l| 
for more details) 



da' 



Til-P)Jo {a -a') 



jr^d^'fia'), (20) 



(a^/)(^) 



and < /3 < 1. Formally, Eq. (|T9)) is obtained by re- 
lating k to the Laplace momentum s of the diffusion 
time a. Concretely, the Laplace transform of an ordinary 
derivative of a function f{a) is given by C[da- f {(j)]{s) = 
s£[/(cr)]-/(0). Applying this to Eq. ((it]) and identifying 



S^ + 2 



(21) 



then modifies the power of s appearing in the right-hand 
side of the Laplace transform and yields 

s-^ {sC[P{x, x', a)] - P(.T, x', 0)} - V^P(x, x' , a) = 0. 

(22) 
Comparing this expression to the definition of the 
Laplace transform of the Caputo derivative. 



C[d"f{a)]{s) = sf'C[f{<j)]-sP-^f{0). 



(23) 



one is led to the RG-improved diffusion equation Eq. (fTOI) 
with /? given by Eq. ((14)) . Notably, the Laplace transform 
also fixes the type of fractional derivative appearing in 
the diffusion equation, since, e.g., a Riemann-Liouville 
fractional derivative would give rise to a different form of 
Eq. ((23]). Also, using the Caputo fractional derivative en- 
sures that the evolution can be formulated as a standard 
Cauchy problem. 

Following [42, |4j] , the general solution of Eq. ((25|) can 
be written in terms of a modified Laplace transform 



Jo 



Pfi{r,a)= / dsAp{s,a)Pi{r,s), 



(24) 



where 



and 



Fi(r,s) 



(47rs)^ 



(25) 



is the sohition of the standard heat-equation (j3]) in d 
dimensions. Constructing A^(s,i) with the help of an 
inverse Laplace transform leads to a modified one-sided 
Levy distribution, which guarantees the existence and 
positivity of P^ (r, a) as long as Pi (r, s) is a proper prob- 
ability density function [43 |. 

Explicitly, for an initial condition in the form of a delta 
distribution the kernel Ap(s,a) can be represented in 
terms of a Fox function, 



A/sis, a) ^ J-^hIs 



,i//3 



(1,1) 
(l,l//3) 



(26) 



Utilizing the series representation of H^ \ , the kernel can 
be written as an infinite sum: 



Afi{s,(T) 



{-ly 



s ^ r(l - /3 - Bn)T{l + n) Kcr^ 



1+n 



(27) 
Substituting this series into Eq. (|24l) and switching to the 
integration variable u ~ s/a^ , a straightforward compu- 
tation shows that the spectral dimension (is(c) resulting 
from Eq. (|19p is again independent of a and given by 
Eq. PT)) . As in the case of nonlinear time, quantum ef- 
fects lead to subdiffusion: 



2d 



/fractional 



r(i + /?) 



(28) 



Comparing this result with Eq. ([T5|) . we observe that the 
difference in the coefficient can be removed by a rescaling 
of the diffusion time cr. Note that the two cases can be 
distinguished by considering either the probability den- 
sity function itself or its higher moments, which carry 
nontrivial information because the process is not Gaus- 
sian. 

We now focus on special cases in which the sum 
Eq. (|27p can be carried out explicitly. For (3 ~ 1/2, one 
obtains 



^1/2(5, cr) 



e •*■' 



while /3 = 1/3 gives 



32/3 
-4i/3(s,ct) = -tttAi 



(ji/3 [(3a)i/3 



(29) 



(30) 



where Ai is the Airy function. Substituting these results 
into Eq. (|24p . the corresponding probability densities are 



Pi/2{r,^) 



as TTT-e *= 



(47rs) 



d/2 



(31) 



Pi/z{r,(j) 



32^ 

rrl/3 



Ai 



ds ■ 



(3<t)1 



(47rs) 



d/2 



(32) 



The first case applies to the UV fixed point regime of 
Asymptotic Safety, whereas the second will play an im- 
portant role for HL spacetimes. For /3 = 1/2, the integral 
in Eq. (|3ip can be performed explicitly and gives rise to a 
linear combination of hypergeometric functions. While, 
for d > 2, the solution P1/2 shows a pole at r — >■ 0, the 
integral of the probability for any finite yet arbitrarily 
small region is finite, thus having a meaningful interpre- 
tation as a probability density. 

The probability density (PT|) actually admits an inter- 
esting interpretation in terms of an iterated Brownian 
motion (IBM) _44 ^]0 This process is defined by a nor- 
mal Brownian motion evaluated "in random time" . More 
precisely, given two independent Brownian motions Pi, 2, 
AriBM(cr) = Pi[|P2(c)|]. For higher-dimensional IBM, the 
process Pi is a d-dimensional vector while P2 remains 
one-dimensional [471. A physical application of IBM is 
diffusion in a crack [49J, which could help to develop an 
intuition for the effective properties of quantum space- 
time. One could picture the effects of quantum gravity 
as turning a smooth background spacetime into an effec- 
tive spacetime where diffusing particles are subjected to 
a randomized diffusion time, and diffuse as if in a "crack" 
formed by quantum fluctuations of the geometry. 

At this stage, it is illustrative to compare the probabil- 
ity densities originating from normal diffusion ([3]), diffu- 
sion in nonlinear time (jisp , and fractional diffusion (jl9[) . 
For the special case d = 4, /? = 1/2 a snapshot taken at 
the diffusion time c = 20 is shown in Fig.Hl This leads us 
to the following central observation: Both anomalous dif- 
fusion processes give rise to the same spectral dimension 
ds(f ). The two different physical processes can, however, 
be distinguished by their diffusion probability. Thus, we 
conclude that a matching spectral dimension does not 
imply that the quantum spacetime found in different ap- 
proaches to quantum gravity is the same. This is a known 
fact from transport theory that should be apprecicated 
in quantum gravity. As we discuss in App. |BJ a variety 
of different stochastic processes, including labyrinthine 
diffusion on fractals, continuous-time random walks, as 
well as fractional and scaled Brownian motion, can ex- 
hibit the same behavior of the mean square displacement 



® This is a special case of a continuous-time random walk (CTRW) 
14811 , where the length of the jump between one site and the next 
has finite variance and the waiting time between one jump and 
the next has a power-law distribution. The diffusion equation 
in this case is given by {da — Vj)P = 0, where the fractional 
derivative is Caputo's. Subdiffusion occurs because the particle 
is trapped in bound states, where it spends far more time than 
in free motion; this type of subdiffusion is then called trapped, or 
due to trapping. 



and, hence, the same value for /3. It is therefore currently 
unclear whether the common result ds = 2 in the UV 
regime of different quantum-gravity approaches is due 
to an accidental degeneracy or, actually, arises due to a 
common physical origin. Matching the spacetime picture 
emerging from different approaches to quantum gravity 
requires the analysis of more refined spectral quantities. 
The probability density of the diffusion process is one 
such quantity, providing a natural generalization of the 
spectral dimension. 

r^ P[r,20] 
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FIG. 4. The probability density P{r,a = 20) in d = 4, 
weighted by r"^ for normal diffusion (red dotted line), power- 
law diffusion [Eq. H15|) with /3 = 1/2, green dashed line] and 
fractional diffusion [Eq. (|19p with (5 = 1/2, thick blue line]. 



IV. DIFFUSION ON ANISOTROPIC 
SPACETIMES 

The strategy of encoding the anomalous diffusion on 
an effective quantum spacetime by a nonlinear diffusion 
time or a fractional time derivative does not lend it- 
self to an immediate generalization to the nonisotropic 
quantum spacetime underlying HL gravity, or any other 
anisotropic setting. In this case, it is an important ob- 
servation [4^, |4y, [50l - [53j that the probability density pip 
can also be obtained as a solution of an ordinary par- 
tial differential equation including higher-order spatial 
derivatives, once suitable source-terms are included. In 
this section we generalize these ideas in order to obtain 
positive-definite solutions for the UV limit of HL gravity 
in three and four spacetime dimensions. The technical 
details of the computation can be found in App. \C\ 



A. HL spacetime with 2 = 2, 3 

In [I4], it was suggested that the diffusion equation de- 
scribing the UV phase of HL gravity is given by 

[d„-d^-{Vly]Pu,ir,t,a)^0, (33) 

where z is the dynamical critical exponent capturing the 
anisotropic scaling in time and space t — > bt,x ^ h^x. 



The value of z is related to the dimension of spacetime 
d = 1 + z. The diffusion implied by Eq. (|55)) is de- 
scribed by a standard diffusion equation in (Euclidean) 
time, while the diffusion on the spatial slices is anoma- 
lous. As explicitly shown in Sec. |TT1 the appearance of 
the higher order spatial derivatives appearing in Eq. p3p 
result in a Pm^ir^t^a) that is not positive semi-definite, 
and thus not a suitable probability density for a diffusion 
process. 

Based on the results of the previous section, we con- 
struct a positive definite Phl('', ijc) as follows. In Ap- 
pendixOit is shown that the positive definite probability 
densities (|3ip and ([5^ also solve the higher-order partial 
differential equation 



[d.~{VlY]Pi,.{r,a)^Sy,{r,a). 



(34) 



The explicit form of the source terms Siu is given in 
Eqs. (jC9p and (JC11|) . respectively. Essentially, these cor- 
respond to delta-function source terms which render the 
diffusion probability positive definite, while leaving the 
spectral dimension unaffcctcdLj However, there is a new 
challenge in interpreting the source term, whose physical 
origin is presently unclear. 

Comparing Eq. (|34p to Eq. (|33p . we sec that wc can 
match the anisotropic diffusion equation by including a 
suitable time-direction. Insisting in the new P^i^{r^t,a) 
being positive definite requires a slight modification of 
the source term, leading to 



[d. - d1 - (ylr\ Pndt.r,<y) = Pi(i,a)5i/,(r,a) . 



(35) 



The solution of Eq. ([35|) is the direct product of the Gaus- 
sian probability density in Euclidean time and the solu- 
tion of anomalous diffusion in space with j3 = 1/ z: 



PuUt,r,a) = Piit,a)P,f,ir,a). 



(36) 



As a product of two manifestly positive semi-definite 
probability densities, Phl is again a positive semi-definite 
quantity. For z = 2, the spatial part of the probability 
density solves the fractional diffusion equation Eq. (|19p 
and the temporal part follows a standard Wiener process 
(see Appendix [U]) . Thus, we can identify the underlying 
stochastic process as iterated Brownian motion in spa- 
tial directions and standard Brownian motion along the 
coordinate temporal direction. 

The spectral dimension implied by the probability den- 
sity ([55)1 is 



D 

1 + — . 
z 



(37) 



^ The higher-derivative equation with 2 = 2 without source-terms 
has been used to describe the diffusion of a test-particle in the 
UV regime of Asymptotically Safe gravity [3,|g|. 



where D denotes the dimension of the spatial shccs. This 
is exactly the behavior expected for the anisotropic space- 
times emerging within HL gravity [l2| . Our construction 
thereby reconciles the previous derivation with a positive 
semi-definite probability density. 

At this stage, it is interesting to pause for the fol- 
lowing observation concerning the relation of HL grav- 
ity and the Monte Carlo simulations carried out within 
CDT. As shown by Eq. (|37)) the anomalous diffusion ef- 
fects observed in HL have their origin in the spatial part 
of the diffusion equation, predicting that the diffusion 
on a spatial slice should lead to d^jspatiai = D/z. For 
3 -f 1 dimensional HL gravity where D = 3, z = 3 this 
implies dg"^] spatial = 1- This, however, seems to be in 
conflict with recent CDT measurements of the spectral 
dimension on spatial slices [2| which, for Z? = 3, reported 
(ig°^ I spatial = 1.5. It would be interesting to study this 
apparent mismatch. 



B. Anomalous diffusion in space and time 

Of course, the technique for generating a diffusion pro- 
cess that is anisotropic in space and time by tensoring 
the corresponding probability densities is not limited to 
the case where the diffusion in time is given by the prob- 
ability density of a standard one-dimensional Brownian 
motion. More generally, one can also consider the case 
where diffusion is anomalous in both space and time with 
the anisotropy captured by two different diffusion coef- 
ficients /3space and /3timo- In this case, a positive semi- 
definite probability density could be obtained from 



Pa, 



tropy(i,r,a) = PA_(t,a) x PA^_(r,a) . (38) 



Adapting the computation outlined in Appendix [Cl it is 
straightforward to construct the corresponding anoma- 
lous diffusion equation satisfied by Panisotropy(i,''7 0'). 
While this situation certainly goes beyond the framework 
of Hofava-Lifshitz gravity where /?time = 1 by construc- 
tion, this case could occur in other theories of quantum 
gravity. Since it is currently unclear if there is a suit- 
able candidate theory of quantum gravity which realizes 
this situation, we refrain from a detailed analysis of this 
possibility at the present stage. 



what is sometimes done in transport theory, i.e., assum- 
ing a time-dependent anomalous exponent. This, for 
instance, was the proposal for multifractional Brownian 
motion [S^, Isil - ISq or other multiscale processes [43 • For 
the remainder of this section we will focus on the mul- 
tiscale geometries emerging within Asymptotically Safe 
Gravity, I, US, [Mi. 

In Q, 5{k) has been determined as a function of 
the scale-dependent Newton's constant and cosmologi- 
cal constant. This function, which can be evaluated nu- 
merically for a given RG trajectory and shows the three 
regimes J = 0,4, 2 as discussed above, is actually a pre- 
diction from the underlying fundamental theory, inde- 
pendent of the precise nature of the fictitious diffusion 
process used to sample the properties of spacctime. Thus, 
the scale-dependence of 5 is actually the same for all dif- 
fusion processes considered in Sec. IIIII This implies, in 
particular, that diffusion in nonlinear time and with a 
fractional diffusion time give rise to the same dg-profile, 
which also agrees with the one constructed in [9[. For a 
typical RG-trajcctory which runs towards positive values 
of the cosmological constant in the IR, the resulting dg (cr) 
is shown as the (blue) solid curve in Fig. [SJ The central 
feature of this curve is its distinguished plateaux where 
ds{a) is approximately constant. From microscopic to 
macroscopic length scales, these are situated at dg = d/2 
(UV regime), dg = '2d/{2 + d) (semiclassical regime), and 
dg = d (classical regime), and can be linked to universal 
features of the gravitational Renormalization Group flow. 
The plateaux are connected by short transition regimes. 
The solution (|T5)) is valid except for these short regimes, 
i.e., in the plateaux S ~ const such that the power-law 
identification holds. 

dsio-) 
-4- 
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V. MULTISCALE DIFFUSION PROCESSES 

So far, we have discussed anomalous diffusion processes 
with 5 ~ const, implying that the spectral dimension 
of the effective quantum spacetime is the same at all 
length scales. In order to connect the quantum space- 
time at sub-Planckian scales, where dg < d, to the classi- 
cal spacetime picture at long distances, where the spec- 
tral dimension should be equal to the topological one 
ds ~ d, S should be promoted to a scale-dependent quan- 
tity. Generalizing S — >■ S{(t) is not much different from 



FIG. 5. Four-dimensional profiles of the scale-dependent 
spectral dimension ds(o-) obtained from Eq. (|17|l with scale- 
dependent S{k) (blue, solid curve) and the multiscale power- 
law diffusion equation p9|) with cti — 0, /3i = 1/2, cr2 = 
10-6, /32 = 1/3, and as = 10"^, /^s = 1 (black, dashed lines) 
Note that there is no dependence on the constants Ci from 
Eq. ^. 



Exploiting the fact that there is presumably no phys- 
ical meaning attached to the transition regimes [9|, this 
plateau structure can be mimicked in different ways. One 



is the simple multiscale power-law diffusion equation 



E^ 



Ji+l 



a)6{a - ffi) 



d 



daP^ 



vl 



P(r,a)=0, 



(39) 

where 9 are Heaviside step functions, Ui denote the tran- 
sition scales between the distinct regimes, and the Ci are 
constants of appropriate dimensionality which are fixed 
by the continuity and normalization of P[r,a). The 
diffusion-time operator could thereby be either the one 
for nonlinear time or fractional diffusion. In both cases 
the solution is given by a sum with step functions which 
glue together probability densities (|16l) or (|24l) for dif- 
ferent Pi. A typical example tailored to the plateau- 
structure obtained from Asymptotic Safety is given by 
the black dashed curve in Fig. [SJ Smoo th p rofiles of 
multiscale geometries can be found in [13, [53 ■ To find 
a smooth Ansatz, the simplest way is to replace a^ in 
Eq. (|15p by a phenomenological profile l'^{<j), so that 



d^ = d 



(9 In cr 



(40) 



The function i contains the hierarchy of fixed scales de- 
termining the onset of the various regimes. 



In practical computations, the average product of the two 
Green's functions is approximated by the product of the 
averaged Green's functions 

((Gf (x,x')G^Jx',.T))) « {{G^{x,x'))){{Gt.{x',x))) . 

(43) 
The Drude-Boltzmann approximation describes the 
probability of the particle to arrive at x' without any 
collisions. 

Moving from quantum mechanics to Asymptotic 
Safety, it is natural to promote the diffusing particle to a 
scalar field living on the quantum spacetime. The Drude- 
Boltzmann approximation constitutes a natural analogue 
of the probe approximation where the diffusing particle 
propagates on a fixed background. Evaluating Eq. (|4ip 
in the coUisionless limit then requires computing the 
Green's function for this field in a situation where fluc- 
tuations of spacetime are integrated out. We expect that 
such a computation can be carried out along the lines 
of Ref . [53 . This should permit a first-principle evalua- 
tion of the probability density without the need of an a 
priori postulate of a diffusion equation, since we expect 
the microscopic theory to determine the type of disorder, 
associated with the quantum fiuctuations of spacetime. 
We hope to come back to this point in the future. 



VI. OUTLOOK: SPECTRAL PROPERTIES 
FROM FIRST PRINCIPLES 



VII. CONCLUSION 



Up to this point, the discussion of the spectral proper- 
ties of the quantum spacetime was based on the quantum- 
improvement of a classical diffusion equation. This con- 
stitutes a valid first step in capturing structural aspects 
of the quantum spacetime. Ultimately, one would like to 
derive the probability density P(r, a) directly from the 
underlying microscopic theory. 

We expect that, within Asymptotic Safety, such a com- 
putation can be carried out by suitably adapting the 
framework of quantum diffusion in a stochastic medium; 
see [53 for a pedagogical introduction. This frame- 
work relates the conditional probability P{x,x',t) for a 
quantum-mechanical particle to propagate from state x 
to x' in time t to the particle's Green's function: 

P(f,f' ,t) =A^ /^^((Gf(x,x')G^_Jx',x))) 



X e 



-[{e+Lo/l^eof /ial + (e-Lo/2^taf /ial] iiot 



(41) 



-,R/A 



Here Ge (x, x') denotes the retarded/ advanced Green's 
function with energy e, a^ is the width of the wave-packet 
in the initial state and the double angular brackets denote 
an averaging over disorder. Assuming that the averaged 
Green's functions depend only slightly on e and that w is 
small compared to CTe, Eq. (|4ip can be approximated by 
the Fourier transformed probability density 



P{x,x' ,uj) 



1 



2ttpo 



{{Gf{x,x')GtJx',x))). (42) 



In this work, we have used the diffusion of a ficti- 
tious probe particle to characterize the spectral prop- 
erties of quantum spacetimes. We have demonstrated 
that many quantum-improved diffusion equations pro- 
posed in the literature actually yield solutions that are 
not positive semi-definite. This deficit is traced back to 
the appearance of higher-derivative operators occurring 
in the "quantum-improved" procedure. Depending on 
the theory under consideration, these operators emerge 
from an RG improvement of the classical diffusion kernel 
(Asymptotic Safety) , the nontrivial scaling of the area op- 
erator in LQG-inspired scenarios, or higher-order spatial 
derivatives capturing an anisotropic scaling of space and 
time (Hofava-Lifshitz gravity). Based on this observa- 
tion, we have proposed new classes of quantum-improved 
diffusion equations, relevant for all theories above, which 
can accommodate the quantum properties of spacetime. 
Their main virtue is to derive the profile of the (scale- 
dependent) spectral dimension ds{cr) obtained in earlier 
calculations while at the same time yielding a positive 
semi-definite probability density. These results place the 
derivation of ds{cr) on solid grou nds while retaining all 
earlier conclusions [a, 0, 0, Hfl, UA, UM ■ This also allows us 
to investigate the probability density as a probe of quan- 
tum spacetime in more detail. Moreover, we explicitly 
identify the underlying stochastic processes of two can- 
didate equations as fractional, scaled, or iterated Brow- 
nian motion. Although the twin problem (i.e., the non- 
uniqueness of the underlying stochastic process to one of 
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our proposals for the diffusion equation) is not solved at 
this stage, the physical insights entailed by the anoma- 
lous diffusion processes could constitute a first step in 
developing an intuitive picture of the underlying quan- 
tum spacetime. 

As wc have shown by the explicit construction of Eqs. 
(fT5|) and p9p . distinct diffusion equations with differ- 
ent underlying stochastic processes (as well as different 
stochastic processes sharing the same diffusion equation) 
can lead to nearly identical profiles ds{(j)- As discussed 
in [iSl ISTI . |60| , multifractional spacetimes can reproduce 
the same profiles, too, even in versions violating ordinary 
Lorentz symmetries. Thus, the spectral dimension con- 
stitutes only a very rough characterization of the quan- 
tum spacetime. A more refined picture has to go be- 
yond the computation of ds{<y)- In this work, we have 
taken the first step in this direction by constructing can- 
didate probability densities for a particle probing the ef- 
fective quantum spacetime of Asymptotically Safe grav- 
ity and Hoi^ava-Lifshitz gravity. We expect that these 
probability densities will, at some point, also be acces- 
sible within other approaches to quantum gravity, fore- 
most the Monte Carlo simulations underlying CDT and 
EDT, thereby offering a novel way to compare the spec- 
tral features of the resulting effective spacetimes. Our 
proposal complements and goes beyond the studies in 
[9|, ll3|, where the comparison was based on dg alone. 

We close with the following cautious remark. Study- 
ing the quantum structure of spacetime based on quan- 
tum improving an ab initio diffusion equation does not 
yield an unambiguous description of the dynamics of the 
probe particle. As we have shown explicitly, different im- 
provement schemes lead to different probability densities. 
Therefore, ultimately a direct investigation of the proba- 
bility density functions derived from first principles will 
be necessary. Based on this data, it can then be checked 
which effective diffusion processes could actually be real- 
ized by a given candidate theory for quantum gravity. 
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Appendix A: Asymptotically Safe gravity-matter 
systems 

For Asymptotically Safe gravity [1^, the scale- 
dependent flowing action Tk[g] implies a scale depen- 
dence of (g^^) on the RG scale k. In the UV fixed-point 
regimqfl 



{g^„ix))k oc k 



(k ^ oo) . 



(Al) 



The identification of fc^ with the eigenvalues of the flat 
Laplacian p^ yields 



[a.-(V2)2]p(x,x',a) = 0. 



(A2) 



The scaling relation (|A1[) can be derived within simple 
truncations of the full effective action, but is actually an 
exact consequence of Asymptotic Safety |7[, i.e., of scale- 
invariance in the UV. The flowing action F^ is given by an 
infinite sum of all operators /„ compatible with the sym- 
metries of the theory with their fc-dependent coupling 
constants: F^ = J2n9n{k) Inif/fj.}^]- Passing to dimen- 
sionless couplings gn{k) = k~'^"gn{k), with the canonical 
dimensionality dn of the coupling, we can reabsorb the 
factor k'^'^ by noting that /^[fc^g^jy] = fc''"/„[5;ii/] under a 
scale transformation of the metric. Now we use the fact 
that in the Asymptotic Safety scenario the dimension- 
less couplings approach a fixed point in the ultraviolet, 
gn{k) -> 5*: Ffc^oo = J2n9* In[k'^9^Li^]- Thus, the /„ only 
depend on the combination k'^g^ui so the solution {g^u)k 
to the effective equations of motion at large k must obey 
a scaling property: If {g^i,y)ka is a solution at feg, then it 
must be related to the solution {gfj,u)k at fc by a simple 
rescaling, {gf,^)k ^ {^1 / k'^) {g f,^) ko ■ Crucially, this does 
not depend on the presence of a particular operator in the 
effective action. Thus, Eq. (jAip is a generic consequence 
of the theory becoming scale-free in the UV. 

Let us include matter degrees of freedom, and consider 
In [gfiL' , ^] , where $ stands for various matter fields (in- 
cluding gauge fields). Using the canonical form of the 
scale transformation of a matter field, the factor k'^" can 
be reabsorbed again. Thus, the argument given in Q 
extends to the more realistic case of Asymptotically Safe 
gravity including arbitrary matter. 

Let us spell out the argument in some more detail: The 
action of gravity coupled to matter can be schematically 
written as Ffc = !]„,,« ff">™(^)^" Iff A"'] ^™[*]- Herein any 
operator can be written as a metric operator multiplied 
by a matter operator. The former can in this case have 
open indices, so the /„ are a wider class of tensor struc- 
tures than in the pure metric case. Again, going over 



Here the notion of a scale-dependent metric makes sense, since, 
using {giiv)k—fOy which is the solution to the full quantum equa- 
tions of motion, one can define the scale k with respect to this 
unique, dynamically determined metric. This allows one to de- 
fine a family of metrics {gfj.v)k^ which describe the effective back- 
ground metric in an effective theory with cutoff scale k. 
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to dimensionless couplings, the effective action becomes 
Tfe = J2n m9n."i{k)k'^'"-'" In[9^„.]J^[<^]■ Here, the di- 



mensionahty of the couphng satisfies dg 



-dr 



V„ 



Accordingly, under a scale transformation of the metric 
ff/if — >■ k'^QiiL' and of the matter fields <& — > fc~''*$, we get 
J™[$] ^ ^-'^■'-J,„[fc-'^*$] and/„[5^,] ^ fc-'*^"/[fc25^,]. 
Therefore, we can again reabsorb the factor k'^s^.-m a^j]^ 
get Tfe = Y.n.,m9n,m{k)In[k^g^,^]J,n[k-'^'^^. In the fixed 
point regime, gn,m{k) — ?► gn,m* and thus the action only 
depends on the combination k^g^u- Again, we arrive at 
the conclusion that 

k"^ 
{9i,u)k = -^{gt,u)ko- (A3) 

In addition, we conclude that, simultaneously, 

W^--fJ(<I')^o- (A4) 

In other words: If we have a regime where all dimen- 
sionful couplings scale according to their canonical di- 
mensionality, then this regime can only be scale-free, if 
the metric expectation value follows the inverse canonical 
Weyl rescaling, as in Eq. (|A3p . 



Appendix B: Degeneracy problems in diffusion 

Here, we will discuss two degeneracy problems from 
transport theory that are important in the context of 
quantum gravity: 

• Different stochastic processes can lead to the same 
result for the spectral dimension and the mean 
square displacement, as can be seen from case (a), 
(b), (c), and (d) below. 

• Different stochastic processes can be described by 
the same diffusion equation and probability density 
function, as is the case for (c) and (d) which are 
relevant for Asympotically Safe gravity. This is also 
known as the twin problem. 

The fact that different models of quantum gravity sport 
the same spectral dimension in the UV has been enticing 
the community into looking for an explanation in terms 
of common properties of the quantum spacetime. As we 
have discussed, the spectral dimension is a very rough 
probe, and the same value can emerge from widely dif- 
ferent settings. However, the degeneracy of the numerical 
value of ds can be resolved when studying the diffusion 
process associated with these effective quantum space- 
times in greater detail. 

A variety of physical, biological, and biophysical sys- 
tems display the phenomenon of subdiffusion (see 3J] for 
references to the literature). In it, the mean squared dis- 
placement (or second moment, or variance) of the test 
particle is, in certain regimes, a power law in diffusion 
time: 



where ^ is a constant and X{(j) (what we called so far the 
stochastic process) is the random variable describing the 
position of the particle at time a. For normal diffusion 
(Brownian motion), /3 = I. 

Anomalous subdiffusion can be explained by trans- 
port models whose details should be determined by ex- 
periments on a casc-by-case basis. The mean squared 
displacement (|Bip can be reproduced by very different 
stochastic processes. We can recognize four main types 
of subdiffusion: 

• (a) Transport on fractals is realized by labyrinthine 
diffusion [3J] : subdiffusion is caused by the geometric and 
topological structure of the fractal and, so to speak, the 
test probe lives in a "crowded" environment and meets a 
number of obstacles and dead ends. Diffusion on a frac- 
tal is, in general, described by a diffusion equation with 
fractional differential operator d^ and an x-dependent 
diffusion coefficient, plus friction (see [38| and references 
therein). 

• (b) In continuous time random walk (CTRW) [48| . 
the length of the jump between one site and the next has 
finite variance and the waiting time between one jump 
and the next has a power-law distribution. Subdiffusion 
occurs because the particle is trapped in bound states, 
where it spends far more time than in free motion; this 
type of subdiffusion is then called trapped, or due to trap- 
ping. In the words of 3J], while the environment of frac- 
tal diffusion is a sort of rough landscape with a number 
of features, the environment of trap models is more akin 
to flat valleys surrounded by high ridges. The diffusion 
equation is of the form {da- — 3^~'^V^)P — 0, where a 
constant diffusion coefficient is replaced by a fractional 
derivative in time [42], or of the form (d^ — V^)P = 0, 
where the first-order diffusion operator is replaced by a 
fractional Caputo derivative of order /3. 

• (c) Fractional Brownian motion (FBM) |33l . l36l 1371 . 
l6l| (see [SJ, [33, [fij] for reviews) and, more generally, dif- 
fusion associated with a generalized Langevin equation 
(GLE) [H m, mi-ill. This class is characterized by a 
complex environment (e.g., a polymer network) of which 
the probe particle is just one part, interacting with the 
rest in a complicated way. These systems generically dis- 
play viscoelastic behavior, and subdiffusion occurs be- 
cause the probe is dragged by the environment. 

Fractional Brownian motion is defined by a GLE with 
a Weyl fractional integral. 



dcrXpsM{(j) = da 

:^da 



' da' ^^-fl^'^ rjia'), (B2) 

r(7) 



where r; is a Gaussian white noise. This process starts 
at past infinity, is Gaussian and its two-point correlation 
function is (Xpbm(ct)Xfbm(^')> = d(kP^-' + k'^T"! - 
\(T — a'\'^^^^), where d is the topological dimension of 
space. For 7 > 1/2, the mean-squared displacement is 



{X^{a)) Qcal^ , 0</3<l, 



(Bl) 



(^pbm('^)) 



2d\a 



27-1 , 



(B3) 
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the Hurst exponent H in the literature is H ~ 7— 1/2. As 
a model of subdiffusion (1/2 < 7 < 1), FBM is said to be 
anti-persistent [6^, i.e., if the particle velocity is positive 
at one given step, it is more probable that it changes 
sign at the next step. When 1 < 7 < 3/2, FBM is 
persistent, and the sign of the velocity at two contiguous 
steps is most probably the same. On the other hand, 
in labyrinthine environments subdiffusion occurs because 
the probe bounces back against obstacles, and it inverts 
its motion. 

FBM is the only stochastic process in ordinary space 
to be self-similar and characterized by stationary incre- 
ments. Self-similarity means that time segments have the 
same behaviour at any time scale, after some rescaling: 



XpBM (Act) = A'^ ' XpBM (cr) . 



(B4) 



Increments are stationary in the sense that their distri- 
bution depends only on the time interval: 

([Xpbm(ct) - Xpbm(<t')]'> = {X^^,,{^ - a')). (B5) 

There exists also a version [33, [S^ with simpler GLE, 
such that Xrl-fbm(o') := [oI^'rjKa) features a fractional 
Riemann-Liouville integral. However, its increments are 
not stationary [33| . 

• (d) Scaled Brownian motion (SBM) [33l435l | . sim- 
ply defined as a Brownian motion with power-law time, 
Xg-Buicr) '■= B{a^). Although there are no known physi- 
cal examples of it, scaled Brownian motion can be used 
as a fitting model of data with anomalous scaling, in 
cases where increments are non-stationary. It also has 
a remarkable interpretation, due to Sokolov [3J], as a 
Gaussian approximation for a "cloud" of continuous-time 
random walks. 

In terms of the white noise rj{a), the associated 
Langevin equation is 



daXs^^i{a) 



a 2 ,y(CT), 



(B6) 



(X|3ja))=2da^ 



clearly different from Eq. (|B2p . SBM is self-similar, 
-'^sBM(Acr) = A'^/^XsBM(cr), and for /3 = 27 - 1 it repro- 
duces both the FBM scaling law (jB4[) and the anomalous 
mean-squared displacement (jB3[) . 

(B7) 

However, contrary to FBM, the increments of SBM are 
not stationary, since their distribution does not depend 
on the time interval only: 

([XsBM(a) - XsBM(a')]') = ([S(a^) - B{a'^)f) 

^{XlJa-a')). (B8) 

This property is crucial to discriminate SBM from FBM 
experimentially, for instance via the first passage time 
distribution [33| . Another difference btween the two pro- 
cesses is that SBM is Markovian, since the scale trans- 
formation cr — > cr'^ preserves time ordering for /? > [33| . 



On the other hand, FBM is non-Markovian, due to the 
memory effect carried by the nonlocal fractional operator 
in its definition. 

The mean squared displacement alone is not sufficient 
to discriminate among these models and one has to look 
also at higher moments. These yield nontrivial informa- 
tion because the statistics of processes of type (a) and (b) 
are not governed by a Gaussian probability density func- 
tion. In general, however, even the diffusion equation it- 
self does not determine the stochastic process univocally, 
and there may exist processes (named twins) sharing not 
only the same variance scaling law, but also the same 
diffusion equation and probability distribution. An ex- 
ample is Eq. (jl5[) . which encodes both a FBM and an 
SBM. In particular, from the diffusion equation alone it 
is not possible to tell whether a process is Markovian or 
not: its soluton P does not carry all the information of 
a stochastic process. 

All known examples in quantum gravity are subdiffu- 
sive, which is why we concentrated on this case. More- 
over, the spectral dimension ds is typically related to the 
exponent /3 by /3 = dg/dn, where dn is the Hausdorff di- 
mension of spacetime. Therefore, what said above can be 
summarized in the two following statements. (A) Quan- 
tum gravity models with the same Hausdorff dimension 
can have the same spectral dimension even if the associ- 
ated diffusion equations are quite different. (B) The twin 
problem may also be present, as in the case of Eq. (|15p . 

Let us remark that the property of being Markovian 
cannot help us presently to distinguish which of the two 
processes is most probably related to Asymptotically Safe 
gravity, or to the description of the continuum limit in 
CDTs. In the first case, one would not expect that, e.g., 
the propagation of a physical particle on the quantum 
spacetime be non-Markovian, but, again, this cannot be 
used as an intuition for our setting. In CDTs, one should 
keep in mind that, although the discrete random walk 
on a CDT background is certainly Markovian, a non- 
Markovianness can still emerge in a nontrivial way in the 
continuum limit. 



Appendix C: Diffusion equations with source terms 

In this appendix we collect the technical details under- 
lying the construction of the partial differential equation 
Eq. ([M)) . The construction manifestly makes use of the 
integral representation of Pjs , Eq. (j24[) , for the two special 
cases /3 = 1/2 and (3 = 1/3. 

The first step in the construction observes that the 
integral representation Eq. (j24p allows one to rewrite 
derivatives with respect to a and r of Pp as derivatives 
with respect to s appearing under the integral. Con- 
cretely, it is straightforward to verify that Eq. (|25p satis- 
fies 



(V^rPi(r,s) = arPi(r,,s) 



n=l,2, 



(CI) 



while the kernels (|29p and (|30l) solve the partial differen- 
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tial equations 



and 



9,Ai/2(s,a)=a2Ai/2(s,a) (C2) 



a„Ai/3(s,a) = -53Ai/3(s,a), (C3) 



respectively. Setting n = 1/(3 = 2, 3, we can use these 
relations to establish that 



daPp=l ds[d^Ap{s,<j)]Pi{r,s) 

i 

dsi-ir[d^Ap{s,a)]Pi{r,s) 



dsAp{s,a)[d^Pi{r,s)]+Sp 



(C4) 



= / dsAfs{s,a)[{WlrP^{r,s)]+Sp 
Jo 

={yirpp+sp. 

The surface terms Sp arising in the partial integration 
are given by 



n-l 



/c=0 

= Y.{-lf^-[d:-'-^A„^{s,amVlfP,{r,s)] 



/c=0 





(C5) 



where Eq. (jCip has been used in the second step. 

We now analyze the surface terms in more detail. The 
first observation is that, because of the exponential fall- 
off of Ap{s, a) for s — > cx), the upper boundary does not 
give a contribution to Sp. On the lower boundary we 
actually have 



lim Pi(r,s) =5^'^\r), 



(C6) 



which is the initial condition for the standard heat equa- 
tion in d-dimensional Euclidean spacetime with diffusion 



time s. Exploiting that Pp(r^ a) is actually subject to the 
same initial conditions, we can rewrite these terms as 



lim Pi(r,s) = P^(r,0) 

s— >-0 



(C7) 



Finally, we have to evaluate the s — )• limit of the kernels 
Ap{s^a) and its derivatives. For the case /3 = 1/2, the 
surface terms contain the zeroth and first derivative of 
Eq. ([29]), given by 



lini Ai/2(s,cr) = -^= 
lim 5sAi/2(s,cr) = 0. 



(C8) 



Substituting these results into Eq. ()C5P then yields the 
source terms for /3 = 1/2: 



51/2 = 



1 



VlPi/2(r,0) 



(C9) 



Here -Pl/2(?', 0) denotes the ct ^ limit of Eq. ((3T|) and 
essentially constitutes a delta-function source term. This 
source term agrees with the one found in [50| . 
Analogously, the case /3 = 1/3 requires evaluating 



limAi/3(s,a)=^^/3^^2/3)' 

lhna.Ai/3(^,^) = - ^2/3r(l/3) ' 
lima2Ai/3(s,a)=0. 

s-i-O 



The resulting surface term is then 
1 



(CIO) 



5l/3 = -T 



ai/3r(2/3) 
1 



(T2/3r(i/3) 



(V^)^Pi/3(r,0) 
V2Pi/3(r,0). 



(Cll) 



Substituting Eq. dCS]) or Eq. ((CTTI) into Eq. dCj]) finally 
yields the partial differential equation p4p with z = 2, 3, 
completing our derivation. 
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